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Taking into account the effect of self-interaction, the dynamics of the quantum fluctuations of the 
inflaton field with \(f> 4 potential is studied in detail. We find that the self interaction efficiently drives 
the initial pure state into a mixed one, which can be understood as a statistical ensemble. Further, 
the expectation value of the squared field operator is found to be converted into the variance of this 
statistical ensemble without giving any significant change in its amplitude. These results verify the 
ansatz of the quantum-to-classical transition that has been assumed in the standard evaluation of 
the amplitude of the primordial fluctuations of the universe. 
PACS number(s): 98.80.Cq, 04.62,+v, 05.40.+j 
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■ The inflationary universe scenarios explain satisfactorily various aspects of the universe, such as the homogeneity, 
isotropy and the amplitude of the primordial fluctuations observed in the microwave background radiation and in the 
f*** \ large scale structure |l| . However, the incompleteness in the discussion of the evaluation of the primordial fluctuations 
^\ . has been pointed out in several papers as is explained below. 

In most of inflationary universe scenarios, the seed of the inhomogeneity of the universe is traced back to the 
quantum fluctuations of the inflaton held generated by the accelerated expansion of the cosmic length scale. In this 
context, the amplification of the quantum fluctuations is characterized by a large amount of squeezing of the state 
vector. A difficulty in interpreting this sate vector exists in the fact that the expectation value of the linear field 
operator does not show inhomogcneities but vanishes while that of the squared field operator becomes very large. 
Usually, one interprets this quantum state as if it were equivalent to a statistical ensemble which has the same amount 
. £^ | of variance that the corresponding quantum operator has in the sense of an expectation value once the scale of the 
fluctuations of our interest exceeds the Hubble horizon scale. Here we refer to the calculation based on this ad-hoc 
classicalization ansatz as "the standard calculation" . 

In order to justify the standard calculation, the stochastic approach to the inflationary universe scenario was 

proposed by Starobinsky |^| and further investigated by many authors ||. The evolution of the order parameter 
that is defined as the expectation value of the spatially averaged field operator was studied. If the physical size of 
the spatial averaging is kept constant, the modes with a large comoving wave number become to contribute to the 
order parameter as the universe expands. The contribution from the newly added small scale modes affects on the 
dynamics of the order parameter as a random noise. As a consequence, the evolution of the order parameter mimics 
the Brownian motion and thus the distribution of the order parameter can be well approximated by a statistical 
ensemble of a classical system influenced by the stochastic noise. 

In these early studies, a free scalar field was investigated and little attention was paid to the effect of interaction. 
However, the fluctuations of a free scalar field in a Friedmann-Robertson- Walker spacetime can be decomposed into 
a set of harmonic oscillators which have time dependent spring constants. Then all of these oscillators are decoupled 
with each other. Hence, no classicalization will be expected in such a decomposed degree of freedom. It follows that 
the each modes, which are added to the order parameter during the inflation, never lose their quantum nature. This 
fact implies that the discussion of the classicalization in the context of the stochastic inflation is incomplete M . We 
think that, in the true theory of the quantum-to-classical transition, the classical nature should be observed even 
if each decomposed mode is considered. In other words, we believe that the change of the number of modes that 
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contribute to the order parameter cannot play an essential role in the appearance of the classical nature of the inflaton 
field. 

In this direction, recently, Lesgourgues, Polarski and Starobinsky [gj proposed a new idea to explain the quantum- 
to-classical transition. They claimed that the equivalence between the large squeezing and the classicality of the state 
can be explained if one discards the tiny contribution due to the decaying mode of the perturbation. However, a free 
scalar field is considered in their approach, and the effect of interaction was not taken into account manifestly. Thus 
we think that their elimination of the decaying mode is still artificial and will not be fully justified although their 
result seems to contain an important suggestion. 

Only recently, the importance of interaction has become emphasized |^-^| ■ As a useful tool to deal with the effect 
of interaction with the environmental degrees of freedom, the closed time path formalism has been developed p0| . 
The total system is divided into "the system" and "the environment" so that the system contains the quantities of 
our interest while the environment does not. By using the closed time path formalism, one can integrate out the 
environmental degrees of freedom to obtain the effective action for the system. Many cosmological issues, such as the 
back reaction to the expansion rate of the universe due to the particle creation, have been investigated m. 

Especially, in the references ||,||, llj , the evolution of the fluctuations of the inflaton filed was examined with the 
aid of the closed time path formalism. It was repeatedly stressed that the correct treatment of the effect of the 
environmental degrees of freedom may relax the problem of the fine tuning of parameters in the inflaton potential. 
Let us consider the A</> model of the chaotic inflation. In the standard calculation, A ~ I0~ 12 is required in order 
to explain the observed value of the primordial density fluctuations JjJ . To the contrary, the authors in the above 
references suggested the possibility that the tuning of the coupling constant can be relaxed to A ~ I0~ 6 . One of the 
main issues of the present paper is to cast a doubt on this statement. 

Here we take a conservative picture of the quantum-to-classical transition based on the paper by Joos and Zch 
|^2| . The basic idea is explained in Sec. 2. There we present two sufficient conditions for the system to possess 
the classical nature. Based on this picture, we investigate the A</> 4 model, which is one of the simplest models of 
the chaotic inflation Q. With the choice of parameter A ~ ICP 12 , the two conditions for the quantum-to-classical 
transition are proved to be satisfied and the fluctuations of the scalar field, 6(j) become of 0(H) as is predicted in 
the standard calculation. Hence, we obtain the conclusion that "the standard calculation" is justified in this simple 
model without any significant modification in the amplitude of the fluctuations as opposed to the prediction given in 
literature [|jT^J|. 

This paper is organized as follows. In Sec. 2 we explain our picture of the quantum-to-classical transition and 
clarify in what situation the system behaves as a classical one. In Sec. 3 we explain our simple model and the basic 
assumptions. Integrating out the environmental degrees of freedom, we derive the effective action for the system. In 
Sec. 4 by using the closed time path formalism, the dynamics of the system is analyzed in a less rigorous but a rather 
intuitive manner. More rigorous treatment based on the quantum master equation for the reduced density matrix is 
provided in Sec. 5. Section 6 is devoted for summary and discussion. 



II. DECOHERENCE MEASURE BETWEEN DIFFERENT WORLDS 

In this section, we explain our picture of the quantum-to-classical transition. At first, let us briefly summarize 
the standard discussion of dccoherence based on the analysis of the reduced density matrix. Suppose that the whole 
system can be divided into two parts: one is the "system", which contains the observables of our interest, and the 
other is the "environment" , which is to be integrated out to obtain the reduced density matrix. We set an initial 
condition for the whole system so that the density matrix is given by that of a pure state and the correlation between 
the system and the environment is absent. Hence the state can be represented as 

p{pi) = l^i/sys ® |^i)cnv cnv(^i| ®sys{£j|- (2-f) 

The reduced density matrix is obtained by taking the partial trace over the environment, 

p(x,x';t) := Tr Gnv aya (x\p(t)\x') sya , (2.2) 

where x and x' are the labels of a complete set of state vectors of the system, say, the eigen state of the 
coordinate. If the interaction between the system and the environment is absent, the loss of coherence docs 
not take place, i.e. p(x,x';t) keeps the form of a pure state, i.e., p(x,x';t) = ay s(x\t) sys sys (t\x'} sys - Hence, 
\p(x, x'; t)\ 2 /p(x, x; t)p{x' , x'; i) = I. However, in the presence of interaction, p(x,x';t) no longer keeps the form 
of a pure state. If \p{x, x 1 ; t)\ 2 /p{x, x; t)p{x' , x'; t) becomes quite small for x ^ x' , we may say that the coherence 
between different states labeled by x and x' disappears. If a; is a continuous parameterization of state vectors, we 
can determine the typical scale Ax such that the coherence between the states labeled by x and x' is lost when 
| a; — x'\ > Ax. As long as we observe the system with a resolution coarser than Ax, we may think that the different 
states have no interference between them and so they can be recognized as independent different worlds. However, in 
the above discussion, the stability of the state through the time evolution was not taken into account. The dynamics 
of the system itself and the effects from the environment cause the broadening of the wave function of the system 
in general. If there is a large amount of broadening, it would be difficult to interpret that the state evolves into a 
statistical ensemble of many different "classical worlds" . 
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Here we propose to take more conservative picture of the quantum-to-classical transition in this paper. As mentioned 
in Introduction, the basic idea is taken from the paper by Joos and Zeh fl^) . We restrict our attention to the case 
in which the interaction term between the system and the environment does not contain the momentum variable of 
the system. We suppose that the sufficient conditions in order for the system to possess the classical nature are the 
following two. 

The first condition is that the total system has a set of wave packets which have a sufficiently peaky probability 
distribution in comparison with the accuracy of our measurement through the whole duration we consider. Say, we 
label the wave packets by their initial peak position, s, as \s;t). Since the interaction term does not contain the 
momentum variable, in an approximate sense each wave packet will be written by the direct product as 

\s;t) = \s;t) sys (g>\s;t) cnv . (2.3) 

Here we should note that the state of the environment |s;i) en v is also labeled by s, since the interaction causes the 
correlation between the system and the environment in the course of their time evolution |l3f| . Strictly speaking, 
what we require here is that not |s; t) but \s; t) sys has a sharp peak and that the peak is stable against the evolution. 
Furthermore, we require that this set of wave packets is complete enough that the initial state, can be decomposed 
into a quantum mechanical superposition of these wave packets, i.e., 

|*) = 5>|*;t>. (2.4) 

s 

This means that each wave packet is sufficiently peaky and stable to be recognized as distinguishable "world" . We 
refer to this condition as "the classicality of the dynamics of the system" . 

The second condition is that the coherence between the different wave packets becomes lost swiftly. We can say 
that this condition is the one that was roughly discussed at the beginning of this section. We refer to this condition 
as "the decoherence between different worlds" . For the total system that satisfies the above mentioned first condition, 
the density matrix at time t will be given by 

p(t) = J2c s c* sl \s;t)(s';t\, (2.5) 

s,s' 

The partial trace over the environment gives the reduced density matrix, 

p{t) =^2c s c*,C(s,s';t)\s;t) sya sys (s';t\ :^V... :/:. (2.6) 

s.s f s.s f 

where we defined the partial reduced density matrix p StS i(t) and C(s, s';t) is given by 

C(s,s';t) := cnv (s;t\s';t) cnv . (2.7) 

If we assume that the initial state is given by a direct product of the state of the system and that of the environment, 
we find C(s, s'; U) — 1. If the interaction between the system and the environment is absent, the loss of coherence does 
not take place, i.e. C(s, s'; t) stays time independent constant for all s and s' . However, in the presence of interaction, 
|s;t)env and \s';t} onv evolve differently and hence C(s,s';t) no longer stay constant for s ^ s', while C(s,s;i) = 1. 
If |C(s, s'; t)\ becomes quite small for s ^ s' , we can say that the diagonalization of p has been occurred. Thus the 
quantity |C(s, s'; t)\ characterizes the degree of decoherence between two different worlds, if the evolution of the wave 
packets is not seriously affected by the environment. 
Here we note that 

\c s c:,C( S , S '-t)\ 2 =Tt(p a ,At)pl,At)) > (2-8) 

In a str ict sense, |s; t) can not be written by a direct product as |s; t) sys £g) |s; t) OILV . So th e definition of C(s, s'; t) given 
in (|jj) is not well defined. However, the expression in the right hand side of Eq. ( |2.S| ) makes sense at any time and 
is expected to give the measure of the decoherence between different worlds. 

As a more convenient alternative of the measure of decoherence, we propose to use 

R = """f M . (2.9) 
max| p s ,s'\ti)\ 

If p 8yS i(t) takes the Gaussian form, R equals to \C(s, s';t)\ besides the small correction due to the determinant factor 
that arises from the Gaussian integral. If R becomes very small for s s' , we recognize that the decoherence between 

different worlds is achieved. 

In some sense, our picture is that of the third quantization of the universe pi or of the decoherence history p5[ . 
Our two requirements for the quantum-to-classical transition may be too strong p6| . But, for the present purpose, 
we do not have to relax these conditions. 
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III. SIMPLE MODEL AND THE EFFECTIVE ACTION 



We consider the following simple model of inflation consisted of a single real scalar field using an approximation 
similar to the one that introduced by Matacz M, although the details are significantly modified. We assume that the 
space time can be approximated by a spatially flat de Sitter space, 



ds 2 = dt 2 - a 2 (t)d 2 x, 



with 



«(*) = 77 



.Hi 



H 



and that the Lagrangian of the inflaton field is given by 



S= / dsa A {s) 



d 3 x 



d$(x, s) 
ds 



2a 2 (s) 



-V(*(X,8)) 



(3.1) 



(3.2) 



(3.3) 



where f2 is a finite comoving volume corresponding to the scale of the fluctuations of our interest and we assumed that 
the effect from outside of this volume can be neglected. For simplicity, we choose SI as a cube with < x, y, z < L. 
The lower boundary of the time integration, ti, is the time at which some appropriate initial condition is set. 
We focus on the dynamics of a spatially averaged field in this volume, 



^ / d 3 x$(x,s). 



(3.4) 



We stress that this averaging is performed only on a finite comoving volume, i.e., a part of the time-constant spatial 
surface. Hence 4>(s) does not represent the homogeneous part of the field <i>(x, s) but it represents the fluctuation of 
scale L. Although it is not essential but, for definiteness, we set the periodic boundary condition on &(x,s). Then 
$(x, s) is decomposed as 



$(x, s) = <p(s) + tp(x, s) 



[<7k cos k • x + </ k sin k ■ x] 



k 



(3.5) 



where k 



2tt 



(i,j,k) is a non-vanishing vector with integer, i,j, and non-negative integer, k. Assuming that the 



potential can be approximated by 



(3.6) 



the action reduces to 



S[(j>, q] 



ds a 3 (s) 



j> 2 ( S )-V(^s)) 



a k 



ds a 3 (s)V" (0( S ))(g£) 2 
k 2 



ds a A (s) 



=: S sys [(j)} + Si n t [<?!>, q] + S c 



m 

r[Q], 



a 2 (s) 



(3.7) 



where dot • means a derivative with respect to t. Expecting that the interaction term is small, i.e., V"(Ms)) is 
so, we consider the perturbative expansion with respect to it. As seen from the notation introduced in Eq. (3.7), the 
spatially averaged field is considered as the system and the other short wave length modes are the environment. 

As was performed in Ref. J9|, we calculate the reduced density matrix for <fi integrating over the environmental 
degrees of freedom, q£. Here we assume that the density matrix of the total system is initially represented by the 
direct product as 



ft = /5(^,^;*i)®f[n/4(9&.9ki / ;*i)- 



(3.8) 
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In order to define the initial quantum state of g£, we suppose that the interaction is switched off before t = tj and we 
set the quantum state of g£ by 

<\0) = 0, (3.9) 
where the annihilation operator a£ is defined by the decomposition of the field operator, 

«k = o£«k + 4X*. (3-10) 
and the positive frequency function, w£, is taken as 

1 e ~ tkn ( i \ 

ul = u := _ u . (I--;- , (3.11) 



2fc a(t) V ^ 

before £ = ti. Here we introduced the conformal time coordinate by r\ := —e~ Ht . Then the reduced density matrix is 
calculated to the second order in Sint [<j>, q] as 



p((p,(f)';t) := I dq p{<j> ,(/>', q,q';t) 

b\ [ Vcj> I V^e X p[i{S sy M-S S y S ^'}}+iSM^<t>'}]p(<t>i,^,t), (3.12) 



and 



iS w [4>,cj>']=-i [ ds A(s)f(s) + i [ ds f ds' A(s)£(s>(s, a') 
t 



where 



and 



^ ds ds' A(s)A(s')v(s,s'), (3.13) 



A( s ) = V"(«£)-F"(0'), S( S ) = i(y"(0) + yV)) 5 (3-14) 



cr k 



Ms , s ') = _ t «( 5 ) 3 a^) a ^^ (K(g)]2[Mo(s , )]2 _ [uo(s)] 2 K(s0]2) ; 

cr k 

= «W 3 ^ / ) 3 ^^ (K(s)] 2 K(s0]2 + K(s)] 2 [uS(s0] 2) . (3 . 15) 
cr k 

In order to evaluate the summation over the discrete modes, we replace it with the integral as 

a k z?r jfe m .„ 

where we take fc m m as a constant of 0(2n/L). Since the proper length scale, a{t)L, of the fluctuations which contribute 
to the formation of the large scale structure or the anisotropics of the cosmic microwave background radiation becomes 
much larger than the Hubble scale, H^ 1 at the end of the inflation era. Therefore for simplicity we concentrate on 
the case in which 

a(t)H = -r?" 1 > femin, (3.17) 

is satisfied. Thus we set the initial condition for the reduced density matrix at a time after the length scale of the 
fluctuations of our interest exceeds the Hubble horizon scale. Or equivalently, we assume that the evolution is free 
from the effect of the environmental degrees of freedom until t = ij. Under this condition, we evaluate the functions 
/(s), /i(s,s') and j/(s,s') approximately. The evaluation of f(s), /i(s, s') and v(s,s') is rather complicated. The 
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details of computation are given in appendix A. As shown there, / and ^-terms contain the ultraviolet divergences 
that require renormalization. After subtraction of these divergences, keeping only the leading terms, we obtain 

/(*) 5T5 lo S ' 



8lT 2 \Pmin(s) 

KS ' S) ' (3 - 18) 

where we introduced p m in{s) '■= a~ 1 (s)k m i n and p ma x(s) '■= a 1 {s)k max . As shown in appendix A, in order to 
subtract the divergent portion in the /x-term of Sif, a integration by part with respect to s' is necessary. Then /i-term 
is evaluated as 

S$ = f dsA(s) f ds' f2 0ii(s,a / )+/io(*,a / ))S(* / )+M6(*,a / )S(* / )l . ( 3 - 19 ) 



where 



a 3 (s)0 / p min (s) \ , 
^( s ' s ) IQ^~ \ H M( s_s )' 



, a 3 (s)MI 

Va{S,S) 167r 2 

. a(s)a 2 (s / )n , ., , ... . 

fi b {s, s') ~ 1 ' \ ' log (k min (n - rf)) . 3.20 

Precisely speaking, the approximated values of fxs are different from the true value by a factor of order unity. Thus 
167T 2 can be replaced by, say, 24n 2 . However, these errors do not change the discussion given below in this paper 
because we just show that the effect of these terms is small and can be neglected. 



IV. AN INTUITIVE INTERPRETATION OF THE EFFECTIVE ACTION 



Here we give a less rigorous but a rather intuitive analysis of the evolution of the averaged field, <f>, in the model 
A$ 4 

with V = -^|— . We defer a more rigorous treatment to the next section. 
The effective action is rewritten in terms of 



4>+ 



(f>A = 4>- <i> ', 



as 



S[<f>,<f>']=n J ds a 3 (s) {0+(«)p A (a) -V'(0(s))¥>a} 

-A f ds <MsVa( s )/(.s) + zA 2 / dscj> + {s)<p A (s) f ds' v(s,s')cP + (s')ipa(s') + S^[<f>,<p'], 



and 



where the cubic or higher terms with respect to ip A arc neglected. 

Since u(s, s') decays fast as s — s' becomes large, here we approximate it as 



where 



E(s) :- 



v(s, s') ~ 5{s - s')a 6 (s)Cl 2 E 2 {s)/(f) 2 + {s) 

Ms) 



a 3 (s)n 



ds' 2v{s,s') 



1/2 



H 3 



72a'7r 2 



(4.1) 



(4.2) 



S<">[^'] = A 2 f ds<f> + (s)ip A (s) f ds 1 [(^( S , S ')+A*a(s,s'))0 2 +(5')+^(s,5')0+( S ')0+(s')l , (4-3) 
Jti Jti 



(4.4) 



(4.5) 



G 



with a' := nk~L. Then the i/-term in the action reduces to 



S u [<j>,<tf] = i\ 2 n 2 / ds a 6 (s)S 2 (s) 



(4.6) 

To manage the effect of the ^-term, we introduce an augsiliary field, £, which represents the Gaussian white noise 

<£(*MO>ea = *(*-*'), (4.7) 
where the subscript EA stands for the ensemble average. Further we introduce the action with noise by 



(4.8) 



S ( [<t>,<f>'\ = fi^ ds a 3 (s) {Ms)va(s) - V'M 8 ))<p A } 

-X j ds ^ + ( S )/( S )-a 3 ( S )r!S( S )e( S ))^A( S ) + ^)[ ( /),0']. 



Then from the fact that 



exp(iS[0,<£']) = (exp (iS«[<M'])> 



EA ■ 



(4.9) 



we can expect that the action with noise, gives the evolution of the field including the effect of the fluctuations 
induced by the environment through the ^-term. Taking the variation of S$ with respect to (p&, we obtain the 

Heisenberg equation for the operator <p+(t). Sandwiching thus obtained Heisenberg equation in between the "bra" 
and "cket" vectors of the initial state, the equation for the expectation value becomes 



A 



ds 



0**(*, «) + Ma(*> *)) (<£+(*)<#(*)) + Hbfa s)(Mt)Ms)^+(s)) = 3(f)£(t). (4.10) 



If we set 



(tp) = follows by construction. Thus we find 

($.(*)) = (<Mi)) 3 + 3(0 + (t))(^ 2 (t)> + (<p(t)} 3 - 
Assuming that the effect of the quantum fluctuations is negligible, i.e., (ip 2 (t)) -C (0+(t)) 2 , we obtain 



(4.H) 
(4.12) 



A 2 



a 3 (t)Q 



i+(*)> / ^ 



(m(t,s) + n a (t,s)) (<Ms)) 2 + 



3(t)£(t). (4.13) 



In the rest of this section, we simply use cf>(t) instead of (4>+(t)). The /- and ^-terms give the correction to the 
evolution of in a deterministic manner while the effect of z^-term gives a stochastic force. 

Now we show that the effect of the /- and /x-terms are negligible under the present condition. The effect of /-term 
is just to modify the potential. Thus we compare the force due to the /-term: 



\f(tm) 

a 3 (t)n 



1 



--^\<j>{t)H z log 

07T 



H 



with that due to the bare potential (the third term in Eq. ( 4.13 )): 



vm)) 



G 



(4.14) 



(4.15) 



Hence, the contribution from /-term can be neglected if <fi 2 H 2 \og(H/p min (t)). Since in the inflationary universe 
scenario the typical value of <j> at the time when the comoving scale of the fluctuations of our interest crosses the 
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Hubble horizon scale during inflation is known to become the order of the Planck scale, m p i. Therefore the above 
inequality holds. Thus we conclude that the /-term can be neglected. 

The effect of /x-term is not a simple change of the potential but has a hereditary one. Again the force coming from 
the \Xi and // a -terms is roughly evaluated as 



A 2 



a 3 (t)n 



ds (n(t, s 
A 2 



/i a (M))0 2 (s) 



< 



16tt 2 
A 2 



<j) 3 (t) I ds' 



i, 



log 



Pmin («5 j 

H 



S(t-s')+H log {kminim ~ff)) 



16tt 



^(t) (HAty 



(4.16) 



where rj t is the conformal time corresponding to the cosmological time t and At is the maximum value that t — ti 
takes. Here we neglected the time-dependence of <p(s) because the slow roll ing c ondition is expected to be sat isfied . 
The details of evaluation of the integral is shown in appendix B. Equation (4.16) is to be compared with Eq. (4.15). 

A 2 

Then we find that the Uj and /i a -terms can be neglected as long as the condition 1 ^> ^ (HAt) holds. This 

I67H 

condition is satisfied for typical values of the model parameters such as A ~ 10~ 12 and HAt ~ 60. 
We turn to the contribution from In the same way, it is evaluated as 



A 



a 3 (t)n 



cj>(t) I ds fib(t, s)4>(s)<j>(s) 



X 2 



< 



16tt 2 
A 2 



4> 2 (t)<t>{t) j ds' ^ log (k min (r)t-r,')) 



327r 



<j} 2 {t)(j){t)At. 



(4.17) 



Also, the details of calculation are provided in ap pendi x B. This term should be compared with the friction term due 
to the cosmic expansion (the second term in Eq. ( 4.13| )): 



Then the ratio of these two terms is evaluated as 



3H(f>(t) 
A 



— (—) 

48tt 2 \H 2 J 



(4.18) 



HAt, and is found to be small. Here we introduced a 



constant, v, as a typical value of the inflaton mass squared: v ~ \<fi(t) 2 /2. For A</> 4 model, typical value for v is given 
by v/H 2 ~ 1/100. Thus we conclude that /Xb-term can be also neglected. Thus we concentrate on the effect of the 
i/-term neglecting /- and /i-terms in the rest of this section. 

Under the condition that the fluctuation, 8<j){t) :— </>(£) — 4>(t)\^=o caused by £(t) is small, the above equation reduces 



6ij>(t) + 3H5<j>(t) + V"{(j>{t))8(j){t) = A5(t)£(t). 
Approximating V"(<fi(t)) and S(t) by constants v and 3, respectively, the equation can be solved as 



ds e~ Al(t - s) £(s) - / ds e" A2(t - s) C(s) 



(4.19) 



(4.20) 



where 



Ai = 



3H + D 



A 2 = 



3H-D 



D = yjm 2 - Av. 



(4.21) 



Since v is much smaller than H 2 for \<p model, we can approximate as Ai ~ 3H and A2 ~ v/3H. Then the fluctuation 
caused by this stochastic field £(t) is evaluated as 



(WO 2 



'EA 



A 2 S 2 / D 2 1 



A 2 



-2\ 2 {t-U) 



2D 2 \3Hv 

\ _ e -2A 2 (t-i 4 ) 



432a 



VK 1 



H 1 



(4.22) 



which means that the u-term broadens the peak width of each wave packet as much as ((<5</>) 2 )ea- Hence, the effect 
is negligible small as long as the width of the packet, (<5</>)wp, is much larger than ((S<P) 2 )ea- In the references 
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IID, this quantity )(S<f))% A ( was interpreted as the real fluctuation which is expected to become classical. Thus our 
interpretation is totally different from theirs. Here we do not further discuss this issue. A rigorous justification of our 
interpretation is provided in the next section. 

The main effect of the ^-term is to reduce the off-diagonal elements of the density matrix and brings the quantum 
state into a decohered one. In the above, we have shown that the effect of the environment is ineffective on the 
evolution of the trajectory of the peak, 4>^,{t), of a wave packet, \I/ ,/,(</>, i), where ip is a label to distinguish different 
wave packets. Here we use the initial peak position as the label of wave packets, i.e., ip = <f>^(ti). 

We set the initial condition at a time, ti, sufficiently after the scale of the fluctuations of our interest, L, crossed the 
Hubble horizon, namely a{ti)L > H . Assuming that the quantum state is not affected much by the environment 
before t = t.j , we take the initial quantum state of the system as a pure state represented by a squeezed vacuum state 
which has a large variance 

(54>) 2 QF ~ H 2 , (4.23) 

which corresponds to a natural vacuum state in the de Sitter space such as the Euclidean vacuum state. We decompose 
the initial wave function into the superposition of the wave packets, t). These wave packets are supposed to 

have a sharp peak at <p^{t) with the width, (<5</>)wp> that is much smaller than (<5</>)qf ~ H. We write the reduced 
density matrix at initial time as 



p[^^;U] = J dip J d^Ci^i^C*^')^^',^. (4.24) 

If the condition, (S<P)q F S> ((S<P) 2 )ea, is not satisfied, this decomposition is of no use because we cannot construct the 
wave packets that do not lose their shape as t ime passes. (Later we find that there is another restriction related with 



the uncertainty relation.) Looking at Eq. (4.6), the evolution of this state under the influence of the environment will 
be approximately given by 

p[4>,<l/;t] ~ J dip J d^'C(iP)^^,t)C*(iP')%(cb',t)exp (-A 2 » 2 a 6 ft)5 2 (V ' '^ ^j , (4.25) 

where we thought of S(i) and (p^{t) — 0^/(£) as constants 3 and ip — ip', respectively. The latter replacement will be 



justified because <j>^, is nearly constant when the slow rolling condition is satisfied. The equation (4.25) indicates that 
the off-diagonal elements are exponentially suppressed when 

(I iM^fXM* 6H 2167T 2 f P min(t) \ 6 fXvY 1 u2 

The factor (Xv/H 2 )" 1 is a large number typically of O(10 14 ), but (p m m(t)/H) 6 becomes extremely small as 
e -6H(t-u) ^ e -360^ Hence, (#</>)q F 3> (<5</ > )d c- Therefore even if we require that the wave packets have a peak that 
is sharp enough to satisfy (<5</>)wp <C (<5</>)q F , it is still possible to choose (50)wp so as to satisfy (<5</>) 2 VP S> (<50)d oc - 
So, if the width of the wave packets, i.e., the coarse graining scale of our view, is appropriately chosen, they lose 
the quantum coherence with each other during the inflation. Thus we conclude that the ^-term leads the quantum 
state of the system effectively into a decohered one, which can be recognized as a statistical ensemble of the states 
represented by wave packets with a sufficiently sharp peak, without any significant distortion of the shape or the peak 
position of each wave packet. 

Before closing this section we must mention the effect related with the uncertainty relation. Here we decomposed 
the initial quantum state of the system into a superposition of wave packets with a small variance with respect to the 
variables in configuration space, <f>. We write the width of the wave packet, (5</>)wp as 1/yT for the later convenience. 
According to the uncertainty principle, the small variance in <fi necessarily indicates the existence of a large variance in 
its conjugate variable, fla 3 (t)(f>(t). This variance may induce a large effect on the succeeding evolution. The possible 
presence of the effect of this kind was first pointed out by Matacz jL7|. The induced variance in (S(f>) 2 will become of 
t ■ 1 2 \ 

J t ds6(j)(s) I . This will be evaluated by using the uncertainty relation, 17a 3 (i)((5(/>(i))uR(^0(i))wp ~ 1, as 

(^)ur := (/ da (SMs)h^j - ^ (TH 2 ) (^0) ■ (4.27) 



O 



Thus we find that this effect is also small compared with ((<5</>) 2 )qf ~ H 2 if £j is set at a time well after the scale of 
our interest exceeds the Hubble horizon scale and unless T is extremely large. Further wc note that 
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(<5^ R rri/ 2 ( P i in {t 



(H) 2 WP I 3 V "'# 3 



(4.28) 



Hence, if ti is set to satisfy Pmin{U) <^ H, we can also choose (<$</>)wp to satisfy (<$</>)wp ^ (<^)ur- The restriction to 
the initial time obtained here is consistent with the general belief that the quantum fluctuations of the inflaton filed 
become classical only after the horizon crossing. 



V. MASTER EQUATION 

In this section we directly study the evolution of the density matrix. We choose one solution of an approximate 
equation of motion of <fi(t) obtained by neglecting the effect of the environment. We denote this classical trajectory 
as 4>(t). Namely, </>(£) satisfies 

${t)+3H<i>(jt)+V'$(jt))=0. (5.1) 

Here we introduce new variables, tp and tp', which represent the deviations of the 4> and <f>' from the classical trajectory 
4> by 

tt> = $ + <p, 4/ = $ + <p', (5.2) 
and assume that tp and ip' are small. Then the effective action S[(j>, 4>'] = iS ays [<f)] — iS sys [(j)'] + iSip[<p, <$'] is reduced to 



S[d>, d>'} = Sl£ ds a 3 ( S )| Q^ 2 ( S ) - ±V"${s))tp\s) + 0(p 3 )^j 



^'\s) - -V"{^)W {s)+0{tp' A ) 
ds A(s)f(s)+i f ds A(s) / ds' i>(s,s')A(s') + SM[<l>, </>']. (5.3) 



The time evolution operator for the density matrix can be obtained by constructing the Hamiltonian recognizing (f> 
and <p' as two different interacting fields. Neglecting the cubic or higher order terms in ip and tp', the Hamiltonian 
corresponding to this action is obtained as 

H(t) = ^^P + (t)P A (t) + na 3 (t)V''m))ip + (t)^ A (t) 

ft 



+/(t)A(f) - A(t) J ds [2 imit, s) + fx a (t, a)) E(s) + fx b (t, s)S(s) + iv(t, s)A(a)J , (5.4) 

where we defined 

p+ := — - — , ip A ■■= tp - tp , (5.5) 



and 



P+ = P + P\ Pa = ^^, (5.6) 



are the conjugate momenta of tp+ and tpA, respectively. We note that P and P' are the conjugate momenta of tp and 
tp' . Since tp and P are quantum mechanical operators, they should be associated with hat," , but in order to keep the 
notational simplicity, we abbreviated them. 

In the above Hamiltonian, there appear Heisenberg operators at a past time. The existence of such operators is 
problematic in solving the evolution of the density matrix. To overcome this difficulty we replace the Heisen berg 
operators at a past time with those at present time by using the solution of lowest order Heisenberg equations p8[ , 
which are given by 

P = -ila 3 V"(^, V= T kP, 

P' = Qa 3 V"^)p' 1 ^--Lp'. (5.7) 

J Jar 
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Approximating V"((j)) by a constant v, we can solve these equations as 

<p+(s) =Ti(8,t)<p+(t)+T 2 (s,t)P±(t), <p A (s) =Ti(s,t)<p±(t) + T 2 (s,t)P+(t), 
and Ti(s,t) and T 2 (s,t), are calculated as 



(5.8) 



T 2 («,t) 



A, 



1 / 



a 3 (t)nD 



e A a (t-a) _ gAi(i-s) 



(5.9) 



where Ai, A2 were constants which were defined in Eq. ( 4.2lj ). For the later purpose, we express ip+(t) in an alternative 
way as 



where the operators 0\{t) and C>2(£) are defined by 



(5.10) 



Oi(t) 



D 
1 



1 



O2W := 5 ( A^+(t) 



a 3 (t)ft 
1 

a 3 (i)Q 



Then /-term, that is the term in H which contains /, becomes 

The /i-term is explicitly written as 

X 2 a 3 (t)fl ( p m in{ty 



(5.11) 



(5.12) 



H U>) 



32tt 2 



■log 



PaW (<K*) 3 + 30(t)V(*)) 



AVa(*)[^(*) + V+(*)] / rfs (/i a (i,s)0 2 (s)+/i 6 (t,s)0(s)0(s) 



A tp&(t)4>(t) / (2/x a (i, s)(/>(s)tp+(s) +/J,b(t, s)(f)(s)ip + (s) + [i b (t, s)4>(s)<p + (s) 



(5.13) 



The first line corresponds to the instantaneous part, //i-term. In the second line, the ratio between the first and 
second term in the round bracket is given by 



(j, a (t,s)(j>(t) 



fl b (t, S )(j)(t) 



a2 0) 4>(t) « 2 ( s ) V v 



(5.14) 



and is found to be much greater than unity. Thus the second term can be neglected. By the same reason, the last 
term in the round bracket in the last line can be also neglected. Then the last line is rewritten by using the relation 
dsTTot ) as 



\ 2 <p&(t)<j)(t)Oi{t) / ds{2^ a {t,s)4>{s)- \ x ^{t,s)4>{s)} 



,Ai(t-«) 



\ 2 VA {t)4>{t)0 2 {t) / ds{2n a (t,s)^(s)-\ 2 fi b (t,s)4(s)}e 



X 2 (t-s) 



(5.15) 



Substituting the explicit form of /j, a (t,s) and //&(£, s), and approximating 0(s) by a constant cj>, we will find that all 
the dominant contribution arises from the terms that contain fx a (t,s). Using the formulas given in appendix B, we 
obtain 



= pLi{i)ip A cj) + n 2 (t)<p A <p + {t) + pL 3 {t)ip A Oi(t) + H4(t)<p A 2 (t), 



(5.16) 



with 



|/il(t)|,|/i 2 (*)|>4(t)| 



< 



A 2 a 3 (t)Q 
16tt 2 ' 



(tf Ai) 2 , 
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< 



16n 2 a 3 {uy 



(HAt) , 



(5.17) 



where again we approximated <f>(t) by a constant <f>. 

The v-teim can be also rewritten by using the relation (|]^) as 



h v = -m(t)(p%{t) - iv 2 (t)(p A (t)P+(t) + o(^ 3 ). 



(5.18) 



The coefficients are evaluated as 



vi(t) ~ X 2 cj){t) I ds v(t,s)<p(s)Ti(s,t) 



\ 2 4> 2 



ft 



a 3 (t)SlH 



3_-A 2 (i-t«) 



1/2 (t) 



A^(£) 

a 6Ht 



a 3 (t)Q X 2 



(is z/(£, s)<fi(s)T 2 (s, £) 



-Aa(t-ti) 



(5.19) 



where we used H{t—U) 3> 1 and \ 2 {t—U) = 0(1). Strictly speaking, the former is not the case for £ ~ f*. However, the 
absolute value of the correct expression does not become much larger than that of this approximate expression. Since 
in the later calculation we will find that the contribution from £ ~ £, does not become significant, this approximation 
is not so bad. 

In the above calculations, we used many crude approximations. But small errors caused by these approximations 
will not significantly affect on the results obtained by the following discussion. 

Putting all the results together, we get the Hamiltonian that does not contain any hereditary terms as 



H 



1 



-PaP+ + tta 3 (t)v(t)(p A ip+ + u(t)tp A 



+A»s(*)¥'a(*)Oi(*) + ^(t)^A(t)0 2 (t) - wi(t)ip A - iv 2 {t)<p A P+, 



(5.20) 



where we defined 



v(t) := V"m)) 



u(t) :=A/(£)^(£) + Ml (£)0(£). 



(A/(£)+ M2 (£)), 



(5.21) 



The master equation can be derived from the above Hamiltonian. In the coordinate representation, it becomes 

d d 

idp/dt = Hp, where P+ and P A in H are to be replaced by — i— and — i— , respectively. We note that the same 

master equation can be derived by mean s o f diff erent meth ods pq |. v (£) is dominated by the first term due to the 
same reason explained around Eqs. (4.14), ( 4.1 5| ) and ( 4.16Q in the preceding section. So we approximate v(t) by a 
constant v as before. Later u(i) is found to be compared with a 3 (t)£lV. Since the former is much smaller than the 
latter because of the same reason, the effect of u(t) is negligible small. However, we keep this term for a while until 
this fact turns out to be manifest. 

In order to solve the evolution of the reduced density matrix p(ip,(p';t), it is better to consider in the k — tp A 
representation, which is defined by 



C(fc, <^a; t) = J dip + exp [-ikf + a 3 (t)n] p(<p + , cp A ; £). 
Then the evolution equation for £ becomes 



dt 



ZHk 



Ok 



d(p A 



Vip A 



dk 



iuip A + pz{t)ip A di + (J-4(t)tp A d 2 - ^if A - a 3 Qiy 2 kip A 



(5.22) 



(5.23) 



where we introduced 



Oi 



d 



d 



a 3 {t)nD 



' dk dip A _ 
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On 



a 3 {t)nD 



Ai 



d 



d 



dk d(fA 



Let us assume the Gaussian form of the density matrix as 



C = A/f exp 



NiXi 



(5.24) 



(5.25) 



where Xi = (fc, (p A ) and A/^ is a time independent normalization constant. Then the evolution equation for the density 
matrix reduces to the following set of equations: 



-M 



dt 



dt 







Ma (M feA - X 2 M kk ) + n± (M feA ~ AiMfefe) 
2^3 {Ma a - A 2 M fcA ) + 2/i 4 (M AA - XiM kA ) 


TO 



where 



M 




N 



N A 



We first consider the evolution of My. For this purpose, we define 

<7i := 3-ff, cr 23 :=3H±D, 



ei 




e 2 



0-3 

of/2 



e 3 



fT 2 

of/2 



Then 



(5.26) 



(5.27) 



(5.28) 



2L> 2 



(e 2 + e 3 - 2ei) 



2D 2 



(2criei - cr 2 e2 - cr 3 e 3 ) . 



(5.29) 



follows. Introducing new parameterization of Mij by 

M = Mid ■ 

the equations for Mj , where j 



M 2 e 2 - 

1,2, or 3, are decoupled like 



dMj 
dt 



= °i M 3 



M 3 e 3 , 



S 



Mj, 



where 





Mi - 2M 3 
-Mi 
2M 3 



fJ>4 



a 3 ilD 



Mi - 2M 2 
2M 2 

-Mj 



(5.30) 



(5.31) 



(5.32) 



Here /x-terms contain Mi but we can solve the above equation as if Smj is a given source term. We solve the above 
equation perturbatively taking A as a small parameter. In this sense, Mj is al so ex panded in powers of A. At the 
lowest order, w e solv e the homogeneous equation without source term in Eq. ( 5.31 ). Then to find the next order 
solution of Eq. ( 5.31 ), we need to solve the equation with the source term, S^ij- At this stage, we can substitute the 



lowest order solution, M_- (t) , into Suj ■ Then Smj can be considered as a given source term. Here we should keep 



in mind the limitation of the present analysis. In deriving Eq. (5.31), only the one loop order correction was taken 
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into account, and the lowest order Heisenberg equation was used to remo ve hereditary terms in the Hamiltonian. 
Thus only the correction up to 0{\ 2 ) is valid. If we baldly solve Eq. (5.31) without regard to this limitation, many 
unphysical pathological features will give arise. 
Formally, the solution is given by 



(5.33) 



where 



SM j (t) = J dse^-^Sujis). 



(5.34) 



The contribution from v-terma can be evaluated without specifying the lowest order solution, M- (t) . The corre- 



sponding inhomogeneous solution is approximately evaluated as 

e -A 2 (t-*i) 



f 8M[ v \t) 
5M?\t) 



9ii V l/(6#) 



(5.35) 



\SAlt\t) 

Next we consider A^-part. In the similar manner, we can get the solution as 



N k (t) 
N A (t) 



1 1 

A 2 Ax 



Ni(t) 
N 2 (t) 



with 



N 3 (t) = e^-^N^ + SNjit), (j = 1,2). 



(5.36) 



(5.37) 



Here SNj (t) is a inhomogeneous solution related to the w-term. In the later discussion, we conclude that the effect of 
these terms can be neglected. Thus for the present purpose, we have only to know their order of magnitude. Hence, 
we roughly approximate the time dependence of u(t) by oc a 3 (t). Then we obtain 



SN 2 (t) 



iu{t) ( -3iT(l 
9H 2 



e-M*-**>)/A a 
1 



(5.38) 



Now we discuss the initial condition for the reduced density matrix. We set the initial condition at a time when 
the size of the fluctuations of our interest, fc f T 1 a(t) = (a~ 1 ri) 1 / 3 a(t), becomes larger than the horizon scale, if -1 , 
where a = (2ir)~ 3 . In our present approximation, the earlier epoch is inaccessible, for we used the evaluation of the 
effective action, i.e., the coefficients /, /i and u, under the assumption that a(t)H 3> k. m i n ~ k c . Here we assume that 
the evolution of the fluctuations for t < U can be approximated by the evolution of a non interacting field. Then the 
wave function will be given by 



^(i^) oc exp 



where 



A- — 



l + i 



A 2 
"2* 



H 

Pc(U 



(5.39) 



(5.40) 



and p c (t) — k c /a(t). 

As mentioned in Sec. 2, we decompose the wave function into a superposition of Gaussian wave packets 



^■4,((f) oc exp 



as 



(5.41) 



(5.42) 



where 
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(5.43) 



As before, ip represents the initial position of the peak of wave packets and it is real. For simplicity, we set T is also 
real. The peak width of the wave packet, T -1 / 2 should be sufficiently small compared with the extension of the wave 
function, yj (<^ 2 )qf ~ H. Here we choose, Y to satisfy 



H 2 T > 



aH 

Pc(U) 



> 1. 



(5.44) 



The latter inequality comes from the fact that the initial condition is set after the scale of the fluctuations of our 
interest becomes larger than the Hubble horizon sale. With this choice of T, we find that F ~ A. 
Further, we assume that the initial condition is set at a time when 



H 3 



« 1, 



is satisfied. This choice of the initial time much simplifies the following analysis. Further we assume 

3 / TT \ 3 



TH 2 < A 



H 



Prnin \pi 



H 



Pmin (^) 



(5.45) 



(5.46) 



We refer to this case specified by Eqs. (5.45) and ( |5.46| ) as case A. In case A, as the restriction to T is mild, we can 
examine the dependence of the evolution of the density matrix on the choice of T. 

Alternatively, instead of the limitation on the initial time Eq. (5.45), we can set a rather strong limitation on the 
initial width of the wave packets as 



- > H 2 T. 

A 



(5.47) 



We refer to this case as case B. 

It will be possible to examine more general cases but the analysis becomes much more complicated. So here we 
restrict our attention to these limited two cases. 



Now the reduced density matrix is also decomposed as Eq.(2.6) 

p(t) oc dip dip' exp 



-iP 2 + —iP' 2 
2 v 2 r 



(5.48) 



where we introduced the partial reduced density matrix, p^i(t), which satisfies the same evolution equation as p(t) 
does and its amplitude describes the coherence between the two worlds labeled by the peaks of the wave packets, ip 
and ip' . The initial condition of the partial density matrix is given by 



p^'ik) oc exp 



{( V -^) 2 + ( V -^') 2 } 



(5.49) 



If we introduce the notation 



Pw(t) = A/" exp 



4e 



m l3 y l y j 



E 



(5.50) 



where = (co + , </?a)- The coefficients m^- and rij are related to those in the corresponding k — A representation, Mj, 
and Ni, as 



TO++ 

n+ = 



Mkk 

ia 3 n Nk 



to + a = ia il— — , toaa = Maa 

TIA — Na — ~TT~~Nk, X - A' C exp 



ML 



n 



+ 



m ++ 



(5.51) 



Thus we obtain 



a 6 (t-)n 2 r 

M kk (U) = y r ; , M kA (U) = 0, M AA (U) = -, 
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N k (U) = ia 3 (U)niP+, N A (U) = ~V>a, 



(5.52) 



where 



/ t±V i i /' 

W+ ■= n » VA ■= W ~ W > 



(5.53) 



Noting that 





we can calculate Mj(ti) as 



01 -1 

-0S/2 1/2 
-0 3 /2 1/2 



Mi(t 
M 2 (ti) 
M 3 (ti) 



A/ fcA 



1 

D 



Ai -1 
-A 2 1 



(5.54) 




(5.55) 



Then the lowest order solution in A is given by M-° (t) — e ai<yt ti ^ l Mj{ti). Roughly speaking, when V is not so large 



and satisfies 



H 3 a 3 (u)n^^»H 2 r, 



(5.56) 



the first term in the right hand side of Eq. ( [5.55 ) dominates. In case B, this is always the case if A ^> p^^U) / H 3 
For the later convenience, we introduce 



r := 



r D 2 a e {u)n 2 



(5.57) 



Then it follws that M 2 {U) ~ a 6 (^)ft 2 /2D 2 f . We note that the order of magnitude of the initial value of the other 
components is the same or smaller than that of M2(U). 

To obtain the next order correction, first we need to evaluate the source term Smj- The order of magnitude of 
z^-term is estimated as 



'Mj 



(5.58) 



while that of /x-term is estimated as 



( 1^1 
c(f) I 

\ D M2\ 
V I ^Ar3l 



< 



X 2 (t, 2 a 6 (t)n 2 H(HAt) 2 
l6ir 2 TH 2 



a 3 (U)/a 3 (t) 



(5.59) 



Now it is clear that the contribution from the /i-term is 4a' (H At) 2 / (f H 2 ) times smaller than that from the z^-term. 
This factor can be set small in both cases A and B by choosing T appropriately. This suppression of the contribution 
from the /z-term is not so trivial. The time dependence of n^(t) is approximately proportional to a 6 (t). Hence, if 

there appears the combination, a~ 3 (t)fi 3 (t)M2(t), in S^li it behaves as a 9 (t) and dominates the source term when 
a 3 (t)/ a 3 (ti ) bec omes exponentially large. The disappearance of this kind of dangerous terms is not manifest in 
Eqs. (^310 and (|||). Here we should note that the contribution of /Lt-term to SM^, is much more suppressed by the 
existence of the factor a 3 (ti)/a 3 (t) compared with that of the z^-term. Thus, SMi and 8M2 might be dominated by 
the /i-term but 5M^ is not. 

Under the conditions for case A or case B, M\{t) and M%{t) are dominated by the inhomogeneous soluti on, SMj(t), 
while the contribution to M2(t) from SM2(t) gives only a small collection to M2(t). Then from Eq. (5.33), we find 



Mi(t) 



\ 2 4> 2 a 3 {t)Q, Hjl-e-W-**)) 
1296^~(I) ' 
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M a (t) ~ e^^-^MaCtj) + 5M 2 (t) ~ "^"V * 3 **-^ + 



4r 



m(t) ~ WJfY (5 ' 60) 

Note that they are all real. 

If we again substitute this solution into S^f] , its time dependence becomes proportional to a 9 (t) while that of 

is given by a 6 (t). Thus it seems that S^r] becomes dominant when a 3 (i)/a 3 (£i) becomes large. However, this does 

not mean the breakdown of the present perturbation scheme. In the present calculation, we used the lowest order 

Heisenberg equation in deriving Eq. (5.31). Thus as mentioned before, such a substitution of SMj into Sj£j is n °t 

allowed. . . . . , . 

Using Eqs. fl5.30| ), ( |5.54| ) and ( pMj ), we finally obtain 



m++W 2M 2 (t) 2M 2 (t 4 ) 

m M (t)= (4D 2 M 2 (t)M a (t) - D 2 M 2 {t)) /2M 2 (t) ~ l8H 2 M 3 (t). 



(5.61) 



It should be mentioned that m+A stays purely imaginary. Thus it does not contribute to the absolute magnitude of 
the density matrix. 

For 7^, with the aid of Eqs. ( |5~5^ ) and (||J), we obtain 



n+(t)~ ia 3 (t)Q 
jia(*)= 



gift) 
2M 2 (t) 

- e Mt-u)^ + + e M*-u) 



*(*) 



3H 2 a 3 (t)Q 



H (l _ e -A 2 (t-tO) 



A 2 



6Ha s (U)n 



2M 2 (t) 



(2M 2 (t)N 2 (t) - 2M 3 (t)N 1 (t) + M x {t){N 2 (t) - N^t))) 



3i? 



^ f ^ 12OTAft(t) 



a 3 (i)ft 



+ e 



x 2 (t-t«). 



2rfMi(t) 

a 3 (i)a 3 (ti)fi : 



(5.62) 



Here we note that in the above evaluations of rriij(t) and rij(t) there is no relevant contribution from <5Mi and <5M2. 
Among SMj a relevant contribution is provided only by SM 3 , which is almost insensitive to the effect of /z-term. 
To see the degree of decoherence, we consider the following ratio 



R 



max|p^/(^)| 



(5.63) 



where max means the maximum value when ip + and ip& are varied. As was discussed in Sec. 2, R is the quantity that 
represents how efficiently the coherence between the two worlds (two wave packets) labeled by, tp and ip', gets lost. 
With this definition, R is evaluated by 



with 



K 



i? = cxp {K(t)~ K(U)), 

_ 1 / jR(n A ) 2 | %(n+) 2 
2 \ mAA w++ 



(5.64) 



(5.65) 



In deriving this formula, we used the fact tha t Aft is a constant and we neglected the small logarithmic correction 
that comes from the Gaussian-integral in Eq. (5.22). Then, after a straight forward calculation, we obtain 



K(ti 



(5.66) 



and 



K(t) 



M 2 (t)M 3 (t) 



r 2 r 



32H 2 a 6 (U)n : 



(5.67) 
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for H(t — ti) 3> 1. In Eq. fl5.67| ), the first and the second terms in the round bracket represent the contribution from 
the corresponding terms in Eq. ( [5.65| ), respectively. It is clear that the second term dominates K(t). 
When the first term in the defining equation of Y dominates, i.e., when Eq. ( |5.56 ) is satisfied, 



< 1. (5.< 



9H 2 a e (U)n 2 9H 2 a 6 (ti)fl 2 

Hence, K{t) is much smaller than K(ti), and R is predominantly determined by K{ti) as 

R~exp(-K{U)). (5.69) 

This means that the coherence between two wave packets labeled by ip and ip' is exponentially suppressed for large 
tp/\ and the typical scale of decoherence is determined by the width of the wave packets. For two wave packets 
with < r _1 , their overlap is large. Hence, it is natural that their coherence is maintained. So the scale of the 
decoherence depends totally on the width of wave packets that we choose. 

In case A, we can set the width of wave packet much smaller so that the second term in the defining equation of T 
dominates. Then 

IT 9H 2 a 6 (ti)Q 2 

K (5.70) 



9H 2 a 6 {U)Q 2 Y 2 
and R becomes 

/ 9H 2 a 6 (U)n 2 2 \ 
R ~ exp ( 4] n rAJ ■ (5-71) 

In this case the typical scale of the decoherence between two different wave packets, (fo/Odcc; wu l be given by 

or 

Thus the decoherence does not occurs if we set the initial wave packets too narrow, i.e, if Y is too large. The best 
choice of Y that minimizes R is given by 

TH> = 3g3fl ^f = 3aH \ (5.73) 

PcV-i) 

With this choice of Y, the typical scale of the decoherence becomes 

2 /„3/ 



4H 2 fp 3 c (t 



This bound mainly comes from the broadening of the wave packet due to the uncertainty rel ation . This fact can be 
understood by seeing that the minimum value of (^V')dcc correspond to the scale given in Eq. (4.27) in the preceding 
section. To see this fact in the present context, we consider the expectation value of, (<^+) 2 , which is expected to 
represent the degree of the broadening of the wave packet. If H 2 Y <?C aH 3 /p^(ti), we have 

((^+) 2 ) ~ 2/m++ ~ 1/f , (5.75) 

and is found to stay almost constant. Instead, if we assume a very large value of T such that violates this condition, 
the broadening of the wave packet at a late time is evaluated as 



<(^ + ) 2 > - f 



YH 2 



(5.76) 



and becomes much larger than the initial width of the wave packet Y 1 . It is easy to see that the minimum broadening 
of the wave packet 
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is also achieved when Eq. (5.73) holds. Comparing this result with Eq. (5.74), we find that the scale of decoherence 
is determined by the broadening of the wave packets. 

However, there is another broadening mechanism due to the i^-tcrm. The fluctuations of the environment behave as 
a stochastic noise, and cause the broa dening of the wave packet. In the last section we estimated the fluctuation in 



responsible for this effect in Eq. (4.22), and we left the task to justify our interpretation of ((5(j)) )ea- Now we return 



to this issue. In evaluating Eq. (5.76 J, we have completely neglected the contribution from the z/-term. As mentioned 
before M%(t) has a small correction, 5M2(t), due to the effect of i/-term. This gives an additional broadening of the 
wave packet; 



2SM 2 (U) A 
a 6 (U)n 2 ~ 216a'7r 2 



H(l 



A, 



(—, 

\H 2 J 



(5.78) 



2 CC , which was evaluated in 



which corresponds to the expression previously derived in Eq. (4.22) with the aid of the augsiliary field and naive 
approximations. At this point, the meaning of the quantity ((<5</>) 2 )ea became transparent. ((<5</>) 2 )ea represents the 
fluctuation caused by the environment. Both in case A and in case B, this broadening effect does not change the width 
of wave packets much compared with the case in which this effect is neglected. This is simply because we restricted 

our attention to the case when M^\t) 3> 5Mi{t) holds for simplicity. 

Then we find that the typical scale of the decoherence is not directly related with {5<f))\ 
Eq. (p~26l) . We find this scale in M AA , which gives exactly the same expression as {8(p)\ cc . To understand the meaning 
of this result, we focus on one diagonal component of the partial reduced density matrix, p^{ip, <p'; t). p^{ip, <p'\ t) also 
has two continuous arguments, ip and ip' . The suppression of the off-diagonal elements of p^(p,p';t) is determined 
by toaAi and p^(p, <p']t) becomes exponentially small for a large pa which satisfies pa 3> {S<j>) dec- Thus we can say 
that the decoherence occurs within a wave packet on the scale of (5(f>)dec- However, for narrower wave packets, the 
classicality of the evolution of the system cannot be maintained through the whole duration of our concern. In the 
present case, this condition of the classicality of the evolution of the system totally determines the minimum width 
of the wave packet. 

Finally we observe the peak location. It is calculated as 



(<?+(*)) 



sR(n+) 



-A 2 (t-«i) 



«(*) 



3H 2 a 3 (t)n 



mi 



-A 2 (*-ti) 



A 2 



(5.79) 



The first term represents the change of the separation of the different trajectories labeled by the initial position of 
the peak, ip+. This is just the term to be attributed to the nature of the model potential. In the present model, 
congruence of the classical trajectories converges gradually. The second term is independent of ip + . This term arises 
because the interaction with the environment was not taken into account when we determine (f>(t). Thus it can be 
interpreted as the correction to <f>(t) due to the effect of the envir onme nt . Th is corr ectio n does not change the motion 
of <j>{t") so much by the same reason that discussed around Eqs. (4.14), (4.15) and (4.16). 



VI. SUMMARY AND DISCUSSION 



We investigated the evolution of the perturbation of the inflaton field with A</> 4 potential after the scale of the 
perturbation exceeds the Hubble horizon scale during the inflation. The effect of the coupling to the smaller scale 
modes through the self interaction was taken into consideration by using the closed time path formalism. That is, the 
smaller scale modes are considered as the environment. The initial condition for the quantum state of the fluctuation 
of the inflaton field was set well after the time of the horizon crossing of the considered mode. The initial state was 
supposed to be given by a pure state density matrix composed of a direct product of the usual Euclidean vacuum state, 
which has the variance of 0(H 2 ). This initial quantum state can be recognized as a quantum mechanical superposition 
of wave packets with a sharp peek. In the present model, we found that the influence of the environment does not 
distort these wave packets much but it extinguishes the coherence between the wave packets with different peek 
positions. The efficiency of the decoherence is so high that the state described by the different wave packets can be 
recognized as completely different worlds. Hence, we can conclude that the initial pure state evolves into a mixed 
state which can be interpreted as a statistical ensemble. 

In the context of the inflationary universe scenario, the primordial fluctuations of the universe are evaluated by 
using an ad-hoc classicalization ansatz such that the expectation value of the squared field operator can be interpreted 
as the amplitude of the variance of the statistical ensemble. In this paper we have shown that this ansatz can be 
verified in a simple model. Thus the result obtained here gives a partial justification of the standard calculation of 
the primordial fluctuations. 

However, the important issue might be whether the inflaton field behaves as classical during the reheating process 
that successively occurs after inflation. This is because the previous study on reheating is mostly based on the 
assumption that the fluctuations have already become classical before the reheating occurs. Thus, it will not be 
necessary that the fluctuations of the inflaton field are kept to be classical during the inflation. In order to prove 



19 



that the standard calculation works, we have only to show that the fluctuations of the inflaton become classical not 
throughout the whole duration of the inflation but at the end of it. In this sense, the condition for the standard 
calculation to be justified might be weaker than the conditions required in the present paper. 

In this paper, we restricted our consideration to a specific model and the effect of the metric perturbation and the 
process of reheating were not taken into account at all. So further study is still required as future work. 
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APPENDIX A: EVALUATION OF F{S), n(S,S') AND v(S,S') 
In this appendix, w e sho w the detai ls of the calculation to obtain the approximate expression for /(s), /i(s, s') and 



/(s, s') given in Eqs. (pls| ) and ( ^20| ) 



We begin with the evaluation of /(s). The explicit expression for /(s) is written down as 



/(*) = 



a(s)£l 
~4^~ 



4 



k 2 rj 2 



(Al) 



This expression is divergent, and needs renormalization. For this purpose, we use the point splitting technique. The 
function, / (s), is basically given by using the Weightman function as f(s) oc J d 3 xG^ (x, s; x, s). Now we regularize 
the expression as J d 3 xG^ (x, s, x', s), where x' is chosen to satisfy \x — x'\ 2 = z 2 . After this regularization, f(s) is 
calculated straightforwardly as 



f(s) 



a 3 (s)fl 

16tt 2 
a 3 (s)Q 

I6n 2 
a 3 (s)n f 1 



/ dp f d(cos8)p (l 



a ip cos 9z 



-a - Pmin(s) + 2h2 - H 2 Ei(-ijw(s)z) + Ei(ip min (s)z) + O(z) 



8tt 2 



^logz-^( 7 -l)- 



Pmin( S ) 



H 2 log(p min (s)) + 0(z) 



(A2) 



where Ei is the exponential integral function. The divergent first term in the last line exists even in the limiting case 
where the background curvature can be neglected. Hence, it should be subtracted by the renormalization procedure. 
The second term is to be attributed to the renormalization of the curvature coupling term, £,R(f> 2 , where R is the 
scalar curvature of the spacetime. Setting an appropriate renormalization condition, we obtain the renormalized /(s) 



a 3 (s)n 
8tt 2 



1 / H 

-^Ln(s) + H 2 log — 



(A3) 



2. n(s,s') 



The function //(s, s') is explicitly written down as 

a 3 (s)a 3 {s')n 



H(s, s') = -i 



4tt 2 
a(s)a(s')n 
l&n 2 



dkk 2 (u*(s) 2 u(s') 2 - u{s) 2 u*(s')< 



2 / ■ \ 2 

I \ I I 



(A4) 



where r\ and v[ is the conformal time corresponding to s and s' . We divide it into two pieces, fJ> r (s, s') and fj, s (s, s'). 
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Hs(s,s') := -i 



a(s)a(s')n 



(e 2lk 



(c.c.) 



(A5) 



is the portion that contains the ultraviolet divergence corresponding to coupling constant renormalization and /x r (s, s') 
is the remaining regular terms defined by /i r (s, s') := /i(s, s') — H s (s, s'). 

The expression for /x r (s, s') is slightly complicated but there is no technical difficulty in its evaluation. After a 
straightforward calculation, we obtain 



fj, r (s,s') 



a(s)a(s')n ( 1 fl 1 



12 \r] rf 



1 

12 



rj r\ 
rj 2 rj' 2 



+ 7 + log [2k min (ri - T)')] 



(A6) 



Here we take into account the fact (3.17) and pick up only the dominant terms in the k m i n — ► limit. However, this 
expression is still complicated. In this paper, we just aim to show that the effect of \x— term is negligible small. For 
this purpose, we can simplify the expression by taking the \q\ <C \r\'\ limit as 



Hr(s,s') 



a 3 (s)CLH 



log (k min (r] -?/)). 



(A7) 



When s ~ s' , this simplified expression is not correct. However, since (x(s,s') vanishes in the coincidence limit, this 
simplification does not underestimate the effect of fi— term. Hence, this simplification will be justified. 

The singular part /x s (s, s') needs regularization as before. Recalling that /i(s, s') is essentially given by the product 
of the Weightman function as 



-i 2 



fi(s,s')x d 3 x G {+) {x,t-0,t') +(c.c), 



(A8) 

we can introduce the point splitting regularization by replacing [G^ (x, t; 0, t')] 2 by G( +S) (x, t; 0, t')G^ (x + e, t; 0, if), 
where |e| 2 = z 2 . Then the regularized expression for /i s (s, s') becomes 



/z s (s, s') = -i 



a(s)a(s')Q 
32^2 



dk 



d(cos 8) ( t 
l ^ 



2ik(ri — r t]' ) — 2ik(rj — rj ) A —ikz cos 9 / a(s) 



(A9) 



In order to subtract the singular term which needs renormalization it is necessary to consider the expression including 
the s' integral 



/ := j ds'T,(s')fi s (s, s') 
Integrating by part, this integral is rewritten as 

I=[a 2 (s')U(s,s')E( S ')] s u - 
where 

e 2ik(r,-r,') + e -2ik{ V -v') 



k min 



— / d{cos9) 
k J -i 



as' 



-ik cos 0z/a(s) 



(A10) 



(All) 



-^L Ms(s , s "). (A12) 



a 2 (s") 



In evaluating U(s, s'), when s =/= s', we can take the z — > limit without any trouble, and easily evaluated as 

U(s, s') ~ -^gj [log (2k mm (v - V')) + 7] • 



(A13) 



On the other hand, when s = s' , U(s, s) becomes 

a(s)fl 



U(s,s) 



16tt 2 



[(1 - 7 - log*) - log(p mi „(s))] , 



(A14) 



which contains the logarithmic divergence corresponding to the coupling constant renormalization. After the renor- 
malization, we obtain 



a(s)fl ( pmin(s) 
f/(s ' s) ^-l6^ l0g 



H 



(A15) 



Combining all the results, finally we get the expression given in Eq. (3.19). 
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3. u(s,s') 



In the same way, the dominant terms in v(s, s') are evaluated as 

aisfais'ffl 



'^P - ^ k 2 dk {u*{sfu{t'f + u(s) 2 u*(t'f) 



a(s)a(s')Q (2 ! + sin [2*^(77-77'); , 



32tt2 \3^r,' 2 kf nm 77-7?' 

The second term becomes proportional to 5(r) — 77') in the k max — ► 00 limit. Since the time dependence of the second 
term is different from that of the first one, we cannot simply discard the second term. However, the important 
quantity for the discussion in the present paper is the integral over s' of the product of is(s, s') and some function 
F(s,s') which is always a smooth function with respect to s — s'. Thus we can conclude that the first term gives the 
dominant contribution, and the second term can be neglected. 



APPENDIX B: APPROXIMATE FORMULAS FOR INTEGRALS 



Here we explain the details of the approximation used in evaluating several integrals. In this appendix, we set 



77 = -e- Hs , 77' = -e- Hs ' 



and rji 



-hu 



First we evaluate the following integral 

h ■■= J ds' H log (k mm (77 - 77')) 

= - I — log {k min (r] - 77 )) . 

Changing the variable into x := (77' — 77)/ 77, the integral becomes 

r(m-v)/v dx r(vi-v)/v dx 



log(-fc min 7y). 



(Bl) 



(B2) 



The second term is explicitly calculated to become H(s — ti)log(—k m i n r]). When — rj)/rj < 1 the first term is 
bounded by 



r(vi-n)/v dx ri dx 

Jo l + x Jo 1 + x 



log a; 



7T 

12' 



To the contrary, when (t]i — r])/r] > 1 the first term is evaluated by 



r(vt-v)/v dx n ^ r(m-v)/v 

Jo l + x 12 A 1 



(vi-v)/v 



dx 

Tx 



log X 



and 



(m-v)/v dx 



Ji 1 



+ x 



\ogx 



< log 



77 



(•»-")/" dx 77,-77 77, 

— — = log log — . 

l + x 77 277 



Thus we find 



Next we consider 



/i| < (if At) 2 . 



(B3) 



(B4) 



(B5) 



(B6) 



ds' H 



a 2 (s) 



log (k mln {r\ - 77')) 



= -rf log (k m in(v -rj[)). 

Jm V 



(B7) 
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After a straightforward calculation by using the same change of variable as before, we obtain 



V i Vi 



log 



Thus wc find 



Finally we consider 



\I 2 \<HAt. 

ds'H^- log (k min (v - V)) 
= r?i / drj'r?' log {k mm (i] - r/)) . 

Jrji 

In the same way, after a straightforward calculation, we obtain 



/a 



7/" 



1 / T) 



3 V»7 



1 log 



»7; 



f 



1 f Vi 



9 V 7 ? 



1 



2\r, 



i 



and it implies 



|J 3 | < -ETAt. 



(B8) 
(B9) 



(BIO) 



1 // V 



3 V»7 



1 log (-kminV) 



(BH) 
(B12) 
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